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ABSTRACT 

Aims. Studying the statistical properties of solar ultraviolet emission lines could provide information about the nature 
of small scale coronal heating. We expand on previous work to investigate these properties. We study whether the 
predicted statistical distribution of ion emission line intensities produced by a specified heating function is affected by 
the isoelectronic sequence to which the ion belongs, as well as the characteristic temperature at which it was formed 
(as found previously). Particular emphasis is placed on the strong resonance lines belonging to the lithium isoelectronic 
sequence. Predictions for emission lines observed by existing space-based UV spectrometers are given. The effects on the 
statistics of a line when observed with a wide-band imaging instrument rather than a spectrometer are also investigated. 
Methods. We use a hydrodynamic model to simulate the UV emission of a loop system heated by nanoflares on small, 
spatially unresolved scales. We select lines emitted at similar temperatures but belonging to different isoelectronic 
groups: Feix and Nevni, Fexu and Mgx, Fexvm, Fexix and Fexxiv. 

Results. Our simulations confirm previous results that almost all lines have an intensity distribution that follows a 
power-law, in a similar way to the heating function. However, only the high temperature lines best preserve the heating 
function's power law index (Fexix being the best ion in the case presented here). The Li isoelectronic lines have different 
statistical properties with respect to the lines from other sequences, due to the extended high temperature tail of their 
contribution functions. However, this is not the case for Fexxiv which may be used as a diagnostic of the coronal 
heating function. We also show that the power-law index of the heating function is effectively preserved when a line is 
observed by a wide-band imaging instrument rather than a spectromenter. 
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1. Introduction 

VUV and X-ray small-scale brightenings are often 
detected in images or spec tra of the solar atmo- 
sphere (e.g. iBerghmans et all [1991 lAletti et all l2000t 



lAschwanden fc Parnelll 120021: iChriste et al.l 12008( 1. These 
observations and theoretical considerations indicate that 
these brightenings may be one manifestation of small-scale 
impulsive heating acting on the unresolved fine scale of the 
corona. Modelling unresolved small-scale brightenings can 
help, for example, to explain the long lifetimes of active 
region loops. Understanding the properties of this heating 
and its contribution to global coronal heating are among 
the most challenging questions in solar physics, important 
starting point is identifying key observables that can diag- 
nose proposed heating functions. 

The study of the frequency distribution of the intensity 
of coronal VUV emission lines is one method for investi- 
gating this problem. The distributions derived from obser- 
vations with a variety of instruments and data sets show 
power laws with in dices mainly between —1.5 and —2 (e.g 
lAschwanden! [2005) . Converting the frequency distribution 
into a distribution of thermal energy for the heating events 
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yields a similar power-law, and this is assumed to represent 
the power-law index of the (unknown) coronal heating func- 
tion. However, behind this method there is the assumption 
that the energy conversion mechanisms do not modify the 
original heating energy distribution, so that the measured 
ev ents distribution is eq uivalent to the heating distribution. 

iParenti et al.l (|2006j) investigated this aspect using a for- 
ward modelling approach: a statistical model of coronal 
heating was posited and properties of the predicted emis- 
sion line intensities and thermal energies were deduced. 
Their main result indicated that only high temperature 
lines, formed when the main cooling mechanism in the loop 
was conduction, are reliable. 

These authors showed that the statistical properties of 
the heating function are conserved by the line intensity sta- 
tistical distribution, only during such a period. The lines 
analyzed by these authors did not belong to the lithium 
isoelectronic sequence. 

As with many other ions formed in the solar atmo- 
sphere, the Li isoelectronic ions are formed principally over 
a narrow temperature interval (w 0.3 dex in log T), al- 
though, unlike most other ions, they also have a tail in 
their temperature distributions that extends to high tem- 
peratures (e.g. Fig. [3]), which leads to the ions sampling a 
wider temperature range than most other ions. 

In the present work, the results of Parenti et al. (2006) 
will be extended to investigate the following questions: 
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— How are the statistical distributions affected if the emis- 
sion lines belong to the lithium-like isoelectronic se- 
quence? 

— How will a statistical distribution measured from a 
wideband imaging instrument compare with a distribu- 
tion from a single spectral line observed by a spectrom- 
eter? 

The work presented here is particularly useful in the 
light of upcoming multi-channel, high resolution solar VUV 
instruments . The Atmospheric Imager Assembly (AIA, 
lGolubll2~006f) is scheduled to be flown on the Solar Dynamics 
Observatory (SDO) in 2009 and has eight distinct chan- 
nels, seven of which are narrow bandpass channels cen- 
tred on strong emission lines. Solar Orbiter is planned 
to be launched in 2015 and VUV imaging and spectro- 
scopic instruments - the Extreme Ultraviolet Imager (EUI) 
and Extreme Ultraviolet Spectomet er (EUS), respectively 
- are part of the str awman pay load (|Hochedez fe al.|[2007t 
lYoung fe al.l l2007a|) . Due to the multiple hot channels 
present on SDO/ AIA, it appears to be a promising in- 
strument, in terms of statistics, for observing small-scale 
heating events. In a more general context, the plasma con- 
ditions in the solar atmosphere can be diagnosed if a set of 
lines formed at each layer of the solar atmosphere, can be 
detected. The coronal Li-like ions, as well as some of the 
Fe ions studied here, have lines at long UV wavelengths 
that can be simultaneously observed with chromospheric 
and transition region lines (e.g. Fexvni and Fexix in the 
SUMER waveband). This is an important consideration for 
future spectrometers where we want access to all layers of 
the solar atmospher e with just two or three spectral bands 
(lYoung fc al.ll2007aft . 

Due to the large number of current and planned VUV 
instruments, the present work concentrates on VUV lines 
and we do not investigate the consequences of the small- 
scale brightening model on X-ray emission. However, to our 
knowledge, soft X-ray microflare studies of this kind have 
not been completed previously, although results from the 
Hinode Soft X-ray Telescope (SXT) may be forthcoming. 

The paper is organised as follows. Section 2 introduces 
the hydrodynamical model used in the simulations. Section 
3 describes how the synthetic spectra are built. Section 4 
provides the results, and conclusions are drawn in Sect. 5. 



2. The model 

The hydrodynamica l model used was described in detail by 
iParenti et al.l (|2006h . and we provide only a brief descrip- 
tion. A coronal loop was modelled by a bundle o f unre- 
solved , identical threads following the approach of ICargiii 
(|l994h . Two thousand threads were used each with a half 
length, L, of 10 9 cm, and a cross-sectional area, A, of 
2 x 10 14 cm 2 . At each time t, a thread was described by 
a single temperature, T, and single density, N. The model 
simulates the heating-cooling cycle of each thread indepen- 
dently. After the thread was heated impulsively, the cooling 
proceeded first by thermal conduction (at high tempera- 
tures) and then by radiation (at low temperatures). During 
the period dominated by conduction, the strand was filled 
with plasma by evaporation from the chromosphere, while 
it drains during the radiation phase. The time evolution 
of the density and temperature in each thread was used to 
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Fig. 1. Probability Distribution Function of the heating 
function used for our simulations. 

calculate the time-dependent synthetic spectra in the entire 
loop. The assumed energy loss was 2.5 x 10~ 4 erg cm" 3 s" 1 . 

The heating function employed by IParenti et al.l (2006) 
wa s derived from a mod el of coronal turbulence developed 
bv iBuchlin et al.l (|2003T ). For the present work, a synthetic 
heating function was used, which simulated a sequence of 
small energy- impulsive events (« 10 21 — 10 24 erg), each 
involving only one thread at a time. The heating function 
had a log- normal distribution with index a = —1.7, over 
the energy range « 10 22 — 10 24 erg (Fig.[TJ. Our use of this 
synthetic heating function is justifiable because power-law 
distributions of observed energy events have been reported, 
and theoretical models have also be en shown to genera te a 
power-law distribution of energies (jParenti et alJ 12006b . A 
synthetic function can also be easily modified, to allow a 
full investigation of model parameters. 

3. The synthetic spectra 

To study the behaviour of the optically thin FUV-UV line 
intensities, it is helpful to consider the column Differential 
Emission Measure (DEM) distribution, which represents 
the amount of material at a given temperature along the 
line of sight h: 



DEM(T) = N' 




where N is electron density along ft at a given temperature 
(T). This quantity is linked to a line intensity (/) for an 
optically thin line by: 

poo 

1=1 A(X) G(T) DEM(T) dT [erg cm" 2 s" 1 sr _1 ] (2) 
J o 

where G(T) is the contribution function that contains the 
atomic physics information for the emission line and is 
predominantly a function of temperature for most allowed 
transitions, and A(X) is the abundance of the element with 
respect to hydrogen. 

In our model, we assumed that the line of sight was 
perpendicular to the loop axes, so that the line of sight 
was given by the diameter of the strand multiplied by 
the number of strands emitting at the given temperature. 
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Under these conditions, the shape of the DEM can be rep- 
resented by the distribution of iV 2 versus T. (For details on 
the dependence of DEM(T) on the model parameters see 
iParenti et al.l (|2006l) and references therein.) 

Figure [5] shows the logarithm of N 2 as a function of the 
logarithm of temperature, integrated over the entire loop 
system (2000 strands) and simulation time (10 5 s). This 
plot helps to interpret our results. 



described by IWilhelm et al.l JT995) . This is one of the 
most useful bands in the VUV because it contains lines 
emitted from the chromosphere to the corona. Another 
SOHO instrument is EIT, which acquires solar images in 
four wavelength bands by usin g multilayer optical coat- 
ings (|Delaboudiniere et al.l ll995). EIS is an ultraviolet spec- 
trometer on board the Hinodc satellite operating in the 
wavelengt h ranges 170-211 A and 246-292 A and was de- 
scribed in lCulhane et al.l (|2007f) . 




As mentioned earlier, in this model the cooling of a 
strand proceeded in two separate phases: the conduction 
phase, where the strand is still at a high temperature; and 
the radiation phase, for lower temperatures. For this rea- 
son, we are able to follow the cooling process in Fig. [21 The 
plasma evaporation from the chromosphere acting during 
the period dominated by thermal conduction fills the high 
temperature part (on the right side from the peak) of the 
N 2 ~T distribution. The lower temperature part of the dis- 
tribution is filled during the radiation phase. The peak cor- 
responds to the temperature location at which part of the 
strands are cooling by condu ction and part by ra diation. 

One of the main results of lParenti et al.l ()2006f ) was that 
the intensity distribution of a spectral line is representative 
of the heating function distribution, if the formation tem- 
perature of the line is within the high temperature region of 
the DEM distribution. This is due to the proximity with the 
moment of heat injection, which produces radiative emis- 
sion that inherits the properties of the heating function. 

For the case presented here, the DEM peaks at about 
log T = 6.6. This implies that lines emitted at lower tem- 
peratures are formed while the cooling in the strand is dom- 
inated by radiation. For log T > 6.9, the cooling is dom- 
inated by conduction. An intermediate condition is found 
between these two temperatures, in which the line emission 
loses all information about the heating function. 

Table Q] lists the lines modelled in the present work, 
and their formation temperatures. These lines are ob- 
servable with the instruments SOHO/SUMER, EIT and 
Hinode/EIS. Where possible, we chose two lines with simi- 
lar formation temperatures but belonging to ions of differ- 
ent isoelectronic sequences. Lines belonging to the Li iso- 
electronic sequence are highlighted. The SUMER instru- 
ment on board SOHO is an ultraviolet spectrometer ob- 
serving in the wavelength range 500-1600 A which was 
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Fig. 3. Contribution functions for the spectroscopic lines 
(top) and the temperature response function for the EIT 
171 and 195 filters (bottom). 



Table 1. List of lines synthesised in this work, and with the 
logarithm of the temperature of their maximum emission. 
These lines can be observed by SOHO/SUMER (Nevm, 
Mgx, Fe xvin, Fexix), SOHO/EIT (Feix, Fexn) and 
Hinode/EIS (Feix, Fexn, Fexxiv) 



log (T max ) 


Li-like 


Others 


5.8 


Nevm 770.41 A 


Feix 171.07 A 


6.1 


Mgx 624.94 A 


Fexn 195.12 A 


6.8 




Fe xvm 974.86 A 


6.9 




Fexix 1118.06 A 


7.24 


Fexxiv 192.03 A 
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Fig. 4. Top: PDFs for the Feix (on the right) and Nevm (on the left) lines. Bottom: PDFs for the Fexn (on the right) 
and MgX (on the left) lines. 



For the theoretical calculation of the line intensi- 
ties, we used th e CHIANTI (v. 5.2, iDere et al.l Il997t 
lLandi et~ 20061) a t omic database and software, adopting 
the Mazzotta et al.l (Il998t) ion fractions and p hotospheric 
element abundances ( Grevesse fc Sauval Il998l) . We simu- 
lated the measured data numbers (DN) from the spectrom- 
eters by assuming observations with a 1 arcsec slit and Is 
cadence and using the standard calibration software avail- 
able in the SolarSoft IDL package. In the case of EIT, we 
simulated the total emission in the channel by integrating 
the emission from all lines in the waveband, as predicted 
by the CHIANTI atomic database, taking into account the 
effective area of the instrument. The Feix 171 channel of 
EIT in particular, is affected by significant contributions 
from the slightly hotter Fex 174.53 A line, which there- 
fore alters its response to the heating function. The 195 
channel has a hot component due to Fexxiv 192.03 A. 



Figure [3] shows the contribution functions G(T) for all lines 
(top) and the EIT response function for the channels (bot- 
tom) used in this work. The G(T) functions of the lithium 
isoelectronic lines in the top plot can be identified by their 
asymmetry towards high temperatures. We anticipate in 
our results that this a high temperature tail will increase 
the number of events with low intensity. 



4. Line intensity distributions 

4.1. "Warm" line properties 

We investigate the intensity distribution of "warm" coronal 
lines, i.e., those lines that form at the average temperature 
of the corona (1-2 MK). 

Figure 0] shows the PDFs for Feix and NeVin 
(log T max = 5.8; top panels), Fexn, and Mgx (log T max 
fi, I : bottom panels). The left panels in the Fig. show the 
PDFs from the Li-like lines. The power-law behaviour of the 
heating function was transmitted to the predicted PDFs of 
each wavelength, although, in each case the index a was far 
smaller than the —1.7 value of the heating function. This 
behav iour is consistent with the findings of iParenti et al.l 
(2006). The minimum value of a was reached for the two Li- 
like lines, because of the increase in the number of strands 
emitting a weak intensity; this was caused by the extension 
of their contribution function to higher temperatures. 

To confirm that high temperature tails produced the 
lower a values for Li-like ions, the contribution functions 
were truncated at high temperature and the calculations 
repeated. The PDFs were found to then have a values that 
were consistent with the Feix and Fexn lines, confirming 
our hypothesis. This comparison demonstrated that the Li- 
like ions had a different response to the heating function 
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than other ions formed at the same temperature, due to 
the high temperature tails of the contribution functions. 

4.2. Imager versus spectrometer 

We then investigated weather the difference in the band- 
width of the instrument could have an effect on the shape 
of the PDF. We simulated the EIT intensities in the 171 
and 195 channels, which were compared with the results 
presented in the previous section (Figure [5]) for the 171 
and 195 emission lines observed by Hinode/EIS. 

Comparing Figs. @] and El the PDFs for the EIS lines 
have indices higher than those derived for the two EIT 
channels. We have therefore demonstrated that observing 
emission lines with a wide band imager can affect the PDF 
of the emission line. We note that the EIS Fe ix line is sit- 
uated at the border of the instrument waveband for which 
the response is low. The predicted DN values are therefore 
very small and the line is not useful in practice for the study 
of small-scale brightenings. 
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Fig. 5. PDF for the counts in the EIT 171 (top) and 195 
(bottom) channels. 



To compare with the spectrometer lines, the intensi- 
ties from the EIT channels were calculated by assuming 
a Is exposure time, which is generally far lower than true 
EIT exposure times. For a more realistic case, we also cal- 
culated the PDFs assuming 60 s exposure times (similar 



exposure times are reached during the high cadence EIT 
"Shutter less" progranQ), and verified that the PDFs re- 
tained the same power-law indices as for the Is case. 



4.3. "Hot" line properties 

By hot lines, we indicate those lines that form at tempera- 
tures higher than the average coronal temperature. 
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Fig. 6. PDFs of Fe XVIII (top) and Fexix (bottom). 

The forbidden lines Fexvin 974.86 A and Fexix 
1118.06 A are important because they appear at long UV 
wavelengths close to cooler lines from the transition region 
and chromosphere. A spectrometer observing these lines 
will therefore observe simultaneously lines from distinct re- 
gions of the solar atmosphere, a feature valuable for under- 
standing the multithermal nature of coronal structures. At 
log T = 6.8 Fe xvni forms, close to the DEM peak (see 
Figs. [H and [3]), where the loop threads can cool by both 
conduction and radiation. Fe xix forms during the cooling 
phase dominated by thermal conduction only (log T = 6.9). 
Figure [5] shows the PDFs for these lines (Fexvin at the top 
and Fe xix at the bottom) . The PDF for Fe XVIII does not 
have a power- law distribution, while Fe xix reproduces well 
the statistical behaviour o f the heating function . The Fe XIX 
line was also modelled bv lParenti et al.l (|2006l ). who found 



See http:/ /side. oma.be/EIT/High-cadence 
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similar results. This finding is then a positive test for ou r 
heating function, which, in contrast to lParenti et al.l ( |2006l) . 
has been completely synthesised. Unfortunately there are 
no strong observed EUV Li-like lines that form at the same 
temperatures as FeXVlll and Fexix. For this reason, we 
cannot directly compare with the results shown in Fig. [Jd] 

As a final test, we compiled the intensity distribution 
for the most prominent line of the Li- like Fexxiv ion, 
which exibits a key flare doublet in the EIS wavelength 
bands at 192.03 A a nd 255.11 A and w hose emission 
peaks at log T = 7.24 (jYoung et al.ll2007bh . If nanoflares 
are responsible for the heating in ARs, we would expect 
faint, non flaring, emiss i on at these temperatures (e.g. 
iPatsourakos fc Klimchu k 2006). In the simulation used 
here, the loop system has low DEM values at these high 
temperatures (Fig. 0), but emission is still found at the 
Fexxiv line, which will be completely formed during the 
conduction phase. For the final reason, we would expect 
that the high temperature wing of the line would not affect 
the slope of the PDF. Figure [7] shows the resulting PDF. As 
expected from the previous results, the distribution follows 
a power-law and the index a is close to that of the heating 
function but in fact a little higher, in contrast to results 
for the other ions that have been studied. We believe that 
this is because the line forms at the very end of the heating 
function energy range, for which the statistics are proba- 
bly too small to be representative (see also the drop of the 
DEM(T) at high temperatures in Fig. EJ). This may be a 
limit to our diagnostic method. 
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Fig. 7. The PDF of Fexxiv. 



5. Conclusions 

In this work, we have modelled the statistical behaviour of 
emission line intensities from a coronal loop, whose heat- 
ing function is distributed as a log-normal distribution in 
energy. We investigated the plasma response for a wide 
range of temperatures, by simulating measurements from 
SOHO/SUMER, Hinode/EIS, and SOHO/EIT lines, an d 
comparing with the earlier work of iParenti et al.l (|2006l) . 
Our results may be summarised as follows: 

— The synthesised heating function used in this work re- 
produces the statistical behaviour of a similar func- 



tion derived by a model used in a previous work 
(|Parenti et al.l l200"rjl) . This would support the use of a 
modified version of this function in further tests with 
heating distributions of different p ower-law indices. 

— We confirm the previous results of IParenti et al.l |2006) 
that the power-law index of the heating function is pre- 
served only by the intensity distributions of hot lines 
(log T > 6.9). 

— The shape of the PDF of the line intensities depends not 
only on the temperature of the line formation but also 
on the iso-electronic sequence of the emitting ion. In 
particular, we have shown that the "warm" Li-like lines 
are inappropriate for this type of diagnostic, due to the 
high temperature tail of their contribution functions. 

— The behaviour of Li-like Fexxiv is not, however, com- 
promised by the high temperature tail of the contribu- 
tion function and the index of the PDF is close to that 
of the heating function. 

— We identify a weakness in the use of imaging instru- 
ments for statistical studies of coronal heating: their 
wavelength bands generally contain more than one 
strong emission line that can affect the PDF. 

In light of the present results, the high temperature 
channels on AIA (Fexxin 133 A) and Solar Orbiter (EUI) 
or the flaring spectroscopic lines on Solar Orbiter/EUS 
could be an important source of information for the coronal 
heating problem. At the same time, sensible contributions 
from other lines in the large band instruments need to be 
carefully investigated. 

The results presented here were obtained using a sim- 
ple, hydrodynamical model, halthough it has been shown 
that the model could reproduce the general properties of 
the plasma well. In the future, we will search for further 
confirmation of our findings by using an evolved version of 
our model. 
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